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Abstract: The equations for Four- Wave-Mixing in a Photonic Crystal 
waveguide are derived accurately. The dispersive nature of slow-light 
enhancement, the impact of Bloch mode reshaping in the nonlinear overlap 
integrals and the tensor nature of the third order polarization are therefore 
taken into account. Numerical calculations reveal substantial differences 
with simpler models, which increase with decreasing group velocity. We 
predict that the gain for a 1.3 mm long, unoptimized GalnP waveguide will 
exceed 10 dB if the pump power exceeds 1 W. 
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1. Introduction 

Slow light (SL) can enable interesting applications in photonics jTI and microwave-photonics 
||2l and it is also expected to enhance nonlinear phenomena. In particular, photonic crystal 
waveguides (PhCWs) present SL propagation 1 3 1 which is predicted to enhance self-phase mod- 
ulation (SPM) miSlIll . One intuitive but powerful picture of this effect represents pulses subject 
to a spatial compression which locally increases the power density |7l, pretty much as cars in 
a highway get closer to each other as their speed is decreased. The dependence of nonlinearity 
on group velocity has been recently observed in PhCW for third harmonic generation |[8l and 
for SPM and three photon absorption ||9l[T0ll. 

Efficient four-wave mixing (FWM) was reported recently in 1.3 mm long, III-V semicon- 
ductor PhCW ifTTl . with a conversion efficiency comparable to that of about 1 m-long highly 
nonlinear fiber lfT2l and about 2 cm-long chalcogenide fiber photonic wire lfT3]| . That result 
confirms the theoretical prediction 11141 of enhanced (FWM) on the basis of the square group 
index scaling factor. However, it must be pointed out that a simple square group index scal- 
ing, to model the FWM enhancement, does not actually take into account several fundamental 
features of the phenomenon. 

The first feature is that the group index in the waveguide is a function of frequency and, 
therefore, it is not the same for the various waves involved in the FWM. Observe in fact in Fig. 
[Tal the numerically calculated ITSll group index (red dots) of a GalnP membrane PhCW. The 
calculation is carried out with pumps placed at = 6 different wavelengths approaching the 
band edge: Xim = [1570 + 10 (M— l)]nm, M — 1, ...,6. The second pump wavelength is X2M = 
Aim + 2nm; the signal and idler frequency are calculated according to ©3^ = CO\m — AoJm, 
('}4M = ft>2M + AOm, where Acom = (Oiiw ~ thus satisfying the FWM frequency condition 
(OiM + (i>2M = (O3M + fi'4M- The waveguide parameters are: a = 480nm (crystal period), d ~ 
0.38fl (hole diameter) and h ~ nOnm (PhC slab height). 

Moreover, the modal superposition of the interacting fields must be also considered. In 
PhCWs this overalap is a function of frequency; in particular, as the frequency approaches 
the bandgap, the mode spreads into the hole region as shown in Fig. [Tal where the intensity dis- 
tribution of the electric field within a cell of the PhCW is shown at three different wavelengths. 
Finally, the PhCWs mode is not constant but periodic along the propagation direction and so, 
differently from other photonic waveguides (slabs, fibers etc.), the careful determination of the 
nonlinear effective coefficients is more complicated U |6l \W\ . 



In this paper, an accurate calculation of all nonlinear effective coefficients necessary to eval- 
uate the FWM interaction in the SL regime of a PhCW, is carried out. The derivation is per- 
formed through a perturbation approach, similarly to fTl\, directly from Maxwell's equations. 
The resulting SPM coefficient corresponds to that obtained by previous derivations IHIJI; the 
XPM coefficient is consistent with the one obtained for multimode propagation llT6l . The FWM 
effective coefficients for a PhCW are determined for the first time to the best of our knowledge. 

2. Derivation of the nonlinear propagation equations 

The starting point of our analysis are Maxwell's equations in the frequency domain where 
the linear permittivity e(r) is a spatial function describing the PhCWs structure and PNL{i^,t) 
accounts for the nonlinear response: 

V X E(r, co) = jcofiH{r, co), V x H(r, co) = -y(oe(r)E(r, co) - jcoPNiir, (o)- (1) 

It is assumed that four signals are propagating in the fundamental TE mode of the PhCW, at 
frequencies that satisfy the FWM condition (Oi + (O2 — (0^ + (O4. The electric and magnetic 
fields are then expanded as: 

J 4 ^4 

E(r,to) = - A;e(r,co;)exp[i(j3,-z)], H(r,co) = - ^ A,-li(r,C(J;)exp[;(j3;z)], (2) 

;=-4,^0 ^i=-4,^0 

where A, are the complex amplitudes and the pairs: e(r,a),)exp(;j3,z),h(r,(a,)exp(ij3,z) are the 
Bloch modes at frequencies ©; (a)_, = —ft),) with j8; = j3((B,) the propagation constant. Bloch 
modes satisfy linear Maxwell's equations individually: 

V X [e,exp(ij3iz)] = ;jLi(0,h,exp(!j3,z), V x [li,exp(ij3,z)] = -!ea),e,exp(;/3,z). (3) 

and also obey the relations: e_, = e*,h_, = — li*. For the sake of simplicity, we have omitted 
the explicit dependence on space and replaced the frequency O; with the subscript /. 

We assume that the nonlinearity is small enough so that complex amplitude. A,- = A, (z), in the 
direction z is slowly varying in comparison to exp(;j3,z) and to the Bloch mode within the cell 
{e,,h,}. Furthermore, the Bloch modes are normalized: /y (e, x +e* x h,) -zdV — 4a. Let us 
stress that, with this choice, |A, p = P, is the active power propagating in the z direction at fre- 
quency (Oi IS). In the following, we will show that this is the natural choice for normalizing the 
Bloch modes when group velocity is substantially different from the phase velocity. Similarly 
as in ref. flTl we consider Eqs. [T] calculated at frequency o,; the second equation is scalarly 
multiplied by e* exp(— ij3,z) and then subtracted from the first, multiplied by h* exp(— ;j3,z). The 
result is integrated over the volume of the PhCW unit cell to obtain: 

^ ^ £ • [e,- X li* + e* X li,] dV = jcot e* exp(-ij3,z) • ^Niir, (Oi)dV (4) 

Here, we also used the hypothesis that A,-, (9A,/iiz are slowly varying functions of z and therefore 
can be taken constant over one unit cell. 

We now introduce the explicit form of the third-order nonlinear polarization PArL(r, o,). Us- 
ing the notation of fTFl it reads: 



P7VL(r, (Oi) = EqX^^^ (r; -co,-, «,-): E;E*E; 



X^^\r; CO,-; cOj, ~0)j, a);):E;E*E,- 



eoz'^' (r; co^iEjE^Ei, (5) 



where i,j,k,l = {1,2,3,4}, with the constraint that {i,j,k,l} are all different in the last term. 
The susceptibility tensor is real because multi-photon absorption can be neglected in GalnP 
at A ~ 1.6jUm. On the right hand side (RHS) of Eq. dS]), the first term represents SPM, the 
summation term XPM, and the last the non-degenerate FWM. Here, third and other harmonic 
generations are neglected by assuming that they will not be phase matched. For the sake of 
brevity the tensor explicit dependence on position and frequencies will also be omitted till the 
end of the derivation. Inserting Eq. (|5]l in Eq.® we obtain: 

dAi iCOiEo 



V 



y \Aj\% [ e*-x^^^-eje*eidV+AjA^Akexp{-icj,APz) U* ■ x^^^'-^j^l^kdV 

JV JV 



(6) 



Here Aj3 = jSs + j34 — j3i — ^2 is the linear phase mismatch and C7, = ±1, where the plus (minus) 
sign applies for / = 3,4 (/ = 1,2). In Eq. |6] the SL enhancement of the nonlinear response is 
hidden in the integrals. In order to make this dependence explicit, we use the identity between 
the electromagnetic energy velocity \e and the group velocit that holds in lossless homoge- 
neous media, in periodic ones |fT9l and in PhCWs 1201 . By projecting the energy velocity along 
the axis unit vector z and using the property that the space-time average magnetic and elec- 
tric energies are equal for Bloch modes II2TI . /io/4/;/h, h* = 1/4/yei •<!*, the following is 
obtained: 

. ^ 1/4 /y(e;xh*+e*xh,0-|Qfy ^ Aa ^ ^ ^ 

""■■^ l/4/^,(£o£.(r)|e,f + Mo|h,f)^/V 2/;,eoe,(r)|e,|2^y e^W, 

Note that energies appearing above are normalized consistently with the choice |A, p = P,-. We 
can now normalize the terms on the RHS of Eq. (|6]l multiplying them by the factors r]f (SPM), 
r\fr\^ (XPM) and r\ir\jr]kr\i (FWM), with 77, = ^2a/(£oW,-Vg,-) = 1, V/. The nonlinear coeffi- 
cients of Eq. |6]are now cast in their canonical form so that the equations governing the FWM 
in PhCW are obtained in a form similar to nonlinear fiber optics Il22l : 

dA- ^ 

—^^iYi\Ai\^A,+2i £ Y,j\Aj\^A,+2iYFiA]AjAte-"'-^l^\ 1,2,3,4, (8) 
with the effective nonlinear coefficients and the relative effective volumes taking the form: 



cV, Vi Wf Jv 3xSL 
n, = ^; - = ^X^erz«(r;co,co,,-a.,,a.,.):e,e;e,.y; (10) 

n2(0ia 1 ^f^kn\^^'^ f £r 



7f/= ^T^; = 11 ITT / —FfT^J ■ 'ir;(Oi;(Oj,-(0,,(Oky.ejeJekdV; (11) 



and where n2 = ^Xxxlx/ {4-£r£Qc) is the bulk, nonlinear refractive index coefficient for a linear 
state of polarization ifTSl . 

So, for each nonlinear effect (SPM, XPM and FWM) we determined: 1) the correct enhance- 
ment factor due to SL; 2) the correct overlap integral. The obtained SPM and XPM coefficients 
are consistent to those previosuly found 11116] |T6l; the FWM coefficient is derived for the first 
time to the best of our knowledge. By observing Eqs. (I9I10I1 11 1 a general rule can be remarked: 
the enhancement factor due to the SL is always given by the geometric mean of the group 
indexes of the waves interacting through the tensor x^^^- 



3. Numerical Results 
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Fig. 1. Comparison with model in ref . 1141 [(a)] Square of the group index and SL en- 
hancement factor S, as a function of the wavelength. In the insets the intensity of elec- 
tric field |e,p of the Bloch mode within an elementary cell is shown at three different 
wavelengths. |(b)| Wavelength dependence of the SL scaling factor and of the effective 
FWM nonlinear coefficient, both normalized to their respective value at An = 1570«ra: 

Ynorm = |7F4 ('^1m)/7F4('^1 1 ) I , S„orm = S{Xim) / S(^\i) ■ 



The theoretical findings of the previous section are applied to the PhCW of Fig. [Ta] To the 
aim of explicitly calculating the tensor products in Eqs. ( 19110111b the theory of ref. [23 1 , that can 
be generally applied to zinc-blend semiconductors (group symmetry 43m), is exploited and then 

xSlx = 2xSly = 2zi^4 = Xxxyy For signal wavelengths 1.52/jot < A < 1.62 fim and given that 
£g ~ 1 .9eV for GalnP, the frequencies (Oi j k.i at which the tensor elements are to be calculated 
satisfy the condition 0.4 < h (Oij^k.i/Eg < 0.43. Then, although Kleinmann symmetry ifTSl is not 

(3) (3) 

satisfied, it is found that the dichroism parameter, defined by Xxxyyl XxxL, can be approximated 



to 0.28 11231 . Thus all tensors can be determined from the above relations, from the knowledge 
of the nonlinear refractive index in GalnP, «2 = 10^'^nz^/W and through ifTsl : 



•Z(3)(r);e,e;e, = i: 



, D Xmnop (r; cor, COj ,-(0i,(0k)e jnelgCkp 

nop 



(12) 



where the summations over the indexes m,n,o,p are made on all possible values of the coordi- 
nate axes {x, y, z} and D is the frequency degeneracy factor lITSl which represents the number of 
distinct permutations of the three frequencies {(Oj, — tO/, 0)^} (D = 3 for SPM, D = 6 for XPM 
and FWM). 

We numerically determined, for all previously defined pump and signal wavelengths the 
Bloch mode electrical field and the dispersion relation flSl, then calculating the FWM SL scal- 



ing factor S = Y[t=\ ^^Jk '■^^ effective nonlinear coefficient 7/74 according to Eq. ( fTTT ). In 
Fig. [Ta]5' can be compared to the square of the group index at the mean frequency, which is 
the approximation used in refs. Ifmri4l : a slight discrepancy appears at the band edge (where 
the SL effect is strong). Note that no particular dispersion engineering of the PhCW has been 
realized on purpose; Aj3 could still be reduced by design, to increase FWM efficiency. To eval- 
uate the effective nonlinear enhancement of PhCW, in Fig. [Tb]the wavelength dependence of S 
is compared to that of 7, revealing that the enhancement of FWM coefficient does not follow 
the pure SL scaling S, the large difference deriving from the decrease in the modal overlap 
integrals. 
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Fig. 2. |(a)| Compaiison of the maximum acliievable parametric gain coefficient, YgynPo (solid 
curves), and the actually achieved gain, g (dashed curves), for two different powers (Pq = 
0.9 W circles, Pq = 1.4W triangles) as a function of SL scaling factor 5. |(b)| Parametric gain 
as a function of Ai for a waveguide L = 1. 3mm long, for three different pump powers. 



By following 1221 we can finally determine the nonlinear phase matching and gain coefficient 
for optical parametric amplification of the wave at frequency 03 : 

2 

Aic = AjS - £ [Yi + 2\i~ jlYj, - 273; - 2Y4i]Pi = AjS + r,„„Po, g = 

where the last terms of each of Eqs. ( fT3] l are obtained for Pi = P2 = Pq/2 and Ygm = YF37f4- 
The effective coefficients Ypm and Ygm describe the strength of the nonlinearity contribution, 
respectively, to the phase matching and to the maximum gain. It is remarkable that, differently 
from fiber optics ll22l . Ypm 7^ Ygm- Figure ( l2al i compares the maximum achievable gain coeffi- 
cient YgmPo to the actual one g, which is limited by the phase mismatch. As the pump power 
increases from Q.9W to I AW the phase mismatch is almost completely canceled by the non- 
linear phase terms and the maximum gain is approached. This fact leads to a dramatic increase 
of the FWM gain G3 = P3{L)/P-i{0) = I + Yg^P^ / g'^ smh^ (gL) which is shown in fig. ^ for 
L = 1.3 mm. 

4. Conclusions 

We have derived the nonlinear equations which describe four-wave mixing in photonic crystal 
waveguides directly from Maxwell's equations. These equations are exact in the limit in which 
we can neglect the changes that the nonlinearity induces in the Bloch modes describing the 
field in the photonic crystal waveguide, a situation by far verified in practice. We demonstrate 
rigorously the explicit dependence of the nonlinear enhancement on the group index and that in 
four-wave mixing (where the fields involved have different group indexes) the dispersive nature 
is rigorously accounted for by the geometric mean of the group indexes of the modes involved. 
Moreover, we demonstrate a substantial correction arising from Bloch mode reshaping in the 
nonlinear field overlap. Finally we account for the tensor nature of the nonlinear polarization. 
As an example, we calculated the gain for a 1.3 mm long waveguide operated at moderately 
small group velocity (v^ > c/20). Particularly, when iig ~ 12 and the coupled pump power 
> IW, the expected gain exceeds 10 dB even if the waveguide is unoptimized. 

The research leading to these results has received funding from the European Community's 
Seventh Framework Programme (FP7/2007-201 1) under grant agreement n. 219299 GOSPEL. 
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